wo 2005/010871 



PCT/AU2004/000M7 



1 



Optical data carrier system 



Field of the Invention 

The present invention xelates broadly to optical data carriers. In particular the present 
invention relates to a method of optimising a sampling function for the holographic storing of 
data in an optical data carrier^ to a method of storing data in a disk-shaped optical data carrier^ 
to an optical data carrier disk and to a method of forming said disk. 

Background of the Invention 

Recent progress in commimication transmission rates has not been matched by a growth 
in data storage capacity. For example 7 years ago the standard system communication speed was 
L6Gb/s for a Tl carrier using copper wire and now it is 320Gb/s (10Gb x 32 channels) using 
optical carriers. This equates to an improvemmt factor in communication transmission of 200. 
In comparison, over the same timeframe the corresponding standard data storage media capacity 
growth has been from approximately 0.8GB for compact discs (CD) to approximately 10GB for 
digital video discs (DVD), representing an improvement factor of only 12.5. This leaves, from a 
cs^acity/speed perspective, the state of development in storage technology lagging behind the 
state of development in transmission technology by a factor of 16. 

Moreover the theoretical storage limit of standard plane surface pixel-based optical 
storage media of approximately 10 Gb per standard CD-sized disk has almost been reached 
using DVD technology. 

Holographic data storage is a term given to technologies in which information is stored 
in form of volumetric structures written inside an optical recording medium. The high bit 
density arising from the three-dimensional volumetric nature of the storage process used in 
holographic data storage may contribute to addressing the shortfall in data storage technologies 
when compared to data transmission technologies as discussed above. 



In at least preferred embodiments, the present invention seeks to provide a novel optical 
data carrier system which is suitable for taking advantage of the three-dimensional volumetric 
nature of an holographic data storage process. 
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Summary of the invention 

In a first aspect the present invmition provides a method of holographically storing data 
as in a series of grating structures including m-level coded elements in an optical data cairier, 
wherein m the method comprising: 

forming a grating sampling function as a direct sum of N partial grating sampling 
functions^ each partial grating sampling function having a phase (^) and amplitude (dn), 
wherein each d„ has m possible values. 

The method can further comprise, conducting an optimisation process to determine a set 
of phases q)^ for which a required maximum refiractive index variation in the optical data carrier 
is related to N^, where 0.5 ^jc^ 1. The required maximum refractive index variation in the 
optical data carrier can be proportional to . Preferably x » 0.5. 

The step of forming a grating sampling function can comprise: 

forming the sampling function as a direct sum of L groups of iV partial grating sampling 
functions, each LxN partial grating sampling function having phases and amplitudes, 
represented by matrices dnu respectively; 

and wherein step of conducting the optimisation process comprises: 

separating flie matrix ^/ into sets of N phases corresponding to the N partial grating 
sampling functions in a given groiip, and one set of Ir phases between the L groups; 

determinuig the sets of phases for each group of N partial grating sampling functions 
from a database having stored therein possible combinations of N coded data elements and 
associated sets of phases; and 

conducting said optimisation process to determine the set of L phases between the L 

groups. 

The optimisation process to determine the set of L phases between the L groups can 
comprises conducting the optimisation process to determine the set of L phases between the L 
groups for which a functional characteristic of the sampling function is miiiinused. 

Preferably the functional characteristic of the sampling function being minimised is a 
mean-square deviation or maximum amplitude. 
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The optimisation process to determine the set of L phases between the L groups, can 
comprise applying a functional analysis to determine the set of L phases between the L groups 
for which a functional characteristic of the sampling function is minimised. The functional 
analysis can comprise a steepest descent (gradient) method. 

The optimisation process to detomine the set of L phases between the L groups may 
comprise approximating the functional characteristic of the sampling function utilising an 
aperiodic autocorrelation function. 

The optimisation process to determine the set of L phases between the L groups may 
further comprise, deriving a gradient of the functional of the sampling function from a 
derivative of the aperiodic autocorrelation function. 

The partial grating sampling functions may comprise one- or multi-dimensional 
functions. 

* « • « * • • 

In a second aspect the present invention provides an optical data carrier configured to 
store data in a plurality of grating stractures, said optical data carrier having at least one data 
reading face through which the grating structures are optically accessible for reading, wherein 
each grating structure comprises a series of m-level coded elements, where m S2, for storage of 
data. 

A required mavin^mri refractive index variation in the optical data carrier is preferably 
related to and wherein 0.5 ^jc^ 1 . Preferably the required maximum refractive index 
variation in the optical data carrier is proportional to N"" and O.S ^ x ^ 1. Preferably a: « 0.5. 

The optical data carrier is preferably disk-shaped. The grating structures may comprise 
one- or multi-dimensional grating structures. 

The optical data carrier can conq)rise a rolled-up material strip in which the plurality of 
grating structures are formed. The optical data carrier may further comprise means for 
maintaining the matmal strip in a rolled-up state. 

The means for maintaining the material strip in a rolled-up state may include, a curable 
material or a mechanical structure. 

In another aspect the present invention provides a method of storing data in an optical 
data carrier, the method comprising the steps of: 
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Storing the data in a material strip, and 

arranging the material strip to form the optical data carrier having a reading fece from 
which the stored data is optically accessible to enable reading the stored data. 

The method may include arranging the material strip to form the data carrier comprises 
spooling the material strip into a disk-shi^ed optical data carrier. 

Preferably the material strip comprises a photosensitive material strip, and the step of 
storing the data comprises inducing refiractive index changes in the photosensitive material strip 
to form grating structures that holographically store the data, wherein a required maximum 
refractive index variation in the grating structures of the optical data carrier is related to N' and 
wherein 0.5 ^ 1 . 

In. a frirther aspect the present invention provides an optical data carrier comprising a 
material strip arranged in a manner such that data stored in the material strip is optically 
accessible from a readmg fece to enable reading of the data stored on the optical data carrier. 
The optical data carrier can be formed by spooling the material strip into a disk. The material 
strip can comprises a plurality of grating structures containing the optical data, wherein each 
grating structure is optically accessible from fbe reading fece. 

The optical data carrier can further comprise means for releasably mai nt ai nin g the 
material strip in Has disk shape. 

In another aspect, the present invention provides a mefliod of forming a disk configured 
to store data in a plurality optical data structures; obtaining a strip-like data carrio: storing 
including the plurality optical data structures; and winding the strip-like data carrier into a disk. 
The Step of obtaining a strip-like data carrier can include, writing the plurality optical data 
structures into a strip-like carrier substrate. The optical data structures are preferably grating 
structures having m-level coded elements where m The method may also include attaching 
adjacent layers of the strip-like data carrier to eadi otiier in the wound disk. 

Brief description of the drawings 

Illustrative embodiments of the present invention will now be described, by way of non- 
limiting example only, with reference to the accompanying drawings, in which: 

Figure la illustrates schonatically a conventional "flaf ' storage of CD/DVD-type disk; 
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Figure lb illustrates schematically a holographic disk embodying the present invention; 

Figure 2 shows a four plots (a) to (d): plots (a) and (b) depict a two-level coding 
sequence in Fourier space; and plots (c) and (d) depict the corresponding amplitude of a grating 
in accordance with an embodiment of the invention without clipping and with clipping 
respectively; 

Figure 3 is a schematic drawing illustrating an optical data carrier writing system 
embodying the present invention; and 

Figure 4 is a schematic drawing illustrating an optical data carrier embodying the present 
invention. 

Detailed description of the embodiments 

Figure lA and IB illustrate schematically the difference between conventional '*£lat" 
storage disks such as CD/DVD-type disks 10 (Figure lA) and an exemplary holographic disk 13 
storage embodying the present invention (Figure IB) . In a CD/DVD-type disk 10 of figure 1 A 
data is stored as a pixels 12 formed on a planar surface 14,. Each pixel, e.g. 12, on the 2D plane 
14 can either reflects a scanning beam to a receiver (bit "1*0 or away from the receiver (bit "0*0. 

In the disc 13 of figure IB, the 2D array of pixels is retained, but instead of a single 
"flat" reflector a volumetric multi-channel grating e.g. 16 is used. The volumetric multi-channel 
grating 16 provides total reflection at certain set of angles indicated by the beams referenced 
with numeral 18), and which correspond to binary "r's and provides total transmission at 
another set of angles referenced with numeral 17, to provide binary "0"s. Alternatively, instead 
of using different measurement angles, a set of measurements at different wavelengths can be 
performed. Thus it can be seen that each pixel e.g. 12, 14 of a CD/DVD-type disk 11 that is 
capable of storing one bit of information can be replaced by a holographic pixel 16 capable of 
storing one bit for each of the predetermined reflection angles or wavelengths. 

Each of the gratings e.g. 16 is directly optically accessible from a reading face 20 of the 
disc 13, i.e. neighbouring gratings, are aligned along the reading face 20 so that no grating 
shadows its neighbour at any of the predetermined group optical reading beam incidence angles 
19. 

It will be shown below that for the standard semiconductor red laser wavelength (650 
mn) such a multi-chamiel grating may have more than 1000 distinct channels meaning that each 
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pixel e.g. 12of the CD/DVD-type disk 10 can be replaced by a holographic pixel 16 storing 
more than 1000 bits of data. 

Estimates related to the simplest (quasi ID) holographic solutions: 

The grating length Az may be estimated bom a simple formula: » 1/ where A^ 
is the spectral width of grating. In this expression, the grating length Azis measured in 
centimetres [cm] and the spectral width AS is measured in special normalised units [1/cm], 
which are related to real life spectral width AS^ (e.g. in [nm]) as A4^ = lOmiQAS^ /A? , where 
Ji is the grating central wavelength in microns. For example for A « 0.65 micron and 
Az - 1 mm we get A^^ - 0.023 nm. 

The maximum number of channels N for a given A and Az is given as 
«0.2A,/A<y^(Az). For example for A «0.65nucron and Az = 1mm we get iNr«6000. 

An estimate for the maximum required refractive index change to create such a grating 
An^ can also be made as follows. 

In ihe first approximation one can express the grating amplitude parameter as 

eo ^^^^^ 

g(z) = -2 yR(S) eix.p(r2iSz)dS , 

where RiS) is reflectivity. Taking 

^R(S) = V^exp(-<5^ /A<5*) 

one may fiirther estimate 

9(z) = -2>/^A<y'v/j?exp(-z*A<5*) or 

This gives an ^ression for An^ of 

A/J^ = IX^n^-^RInl Az . 

For /? = 0.5 , = 0.65 microns, = 1.5 , and Az = 0.1 cm we get A/i^ « 10"' , which 
is a realistic requirement for UV interferometric grating writing techniques. 
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However, it should be noted that this estimate is for a single channel grating only. If an 
N channel grating is to be written the necessary refractive index variation will increase by a 
factor of N, meaning that for any significant N the required maximum refractive index variation 
becomes prohibitively large, e.g. for N=1000, An^ « 1 

Thus if an N-channel grating is to be written an optimisation method is needed to reduce 
the required maximum refractive index variation to within physically acceptable limits. In the 
preferred embodiment of the present invention a dephasing optimisation scheme is used that 
reduces the factor of N by a factor of ^^N . Thus, for // = 1000 , the illustrative embodiment 
presented below gives an extra factor of about 30 or An^ »3-10"^, which corresponds to 
state-of-the-art materials. 

Moreover photosensitive materials capable of producing the refractive index change of 
Aw^ « 10"' or more (together with other general requirements like low shrinkage, good 
mechanical properties, etc) may become available in the ftiture. 

The example embodiments to be described below provide a Fourier data encoding 
approach applicable to higjh cs^acity optical data storage and other fields, which can be 
implemented with readily available materials. 

In the example embodiment, a one-level (identical signals) dephasing approach used for 
multi-chaimel fiber Bragg grating optimisation is modified to accommodate coded signals. 
Instead of using time-consuming fiilly nimaerical algorithms, the problem is simplified by 
calculating the so-called reduced mean deviation Amotion and its derivatives with respect to 
phase values analytically, drastically reducing the amount of numerical calculation. 

A block (convolution) optimization approach is applied, where the total number of 
encoded information elements is fru:torized and optimized in distinct steps. The preset 
inventors have found that, by sacrificing the optimisation quality by a small amount (a few 
percent) substantial speed gains can be made in the numerical calculations. 

As a result, the illustrative embodiment advantageously provides a fast and efficient m- 
level Fourier encoding optimization £^proach. Below it will be described for m=2, that is 2- 
level coded or binary elements, but it will be s^preciated by a person skilled in the art that the 
present invention can be expanded or generalised to higiher-level coded elements fi.e. m > 
2^without departing from the spirit or scope of the present invention. 
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Similar to a single-channel fiber gragg (^rating (FBG), one bit of infonnation can be 
written into a body of photosensitive material e.g- a holographic data disk, by irradiating the 
photosensitive material of the disk with laser light, the intensity profile of which is proportional 
to the grating design, and can be represented as follows: 

qiz) = Kiz)eKp[i(KoZ -h ^(z))] , 

where z is the depth direction of the disk; 
/c(z) is the grating amplitude 
0{z) is the grating phase, and 

Ko is the grating wave number. 
The grating amplitude <(z) is measured in cm'^ and is related to the disk effective 
refractive index modulation An as: . 

ic(z) = ;zAn(z)/(2Aono) , 

where Ao is the grating period and is the average refractive index of the disk. 
Recording LxN bits (iV is the number of bits in a codeword, L is the number of codewords) 
into the disk can be accomplished via periodical modulation (i.e. sampling) of the one-bit 
design: 

qiz) = K(z)exp[iiK^z + 0)]Siz) , 

where Siz) is a periodical sampling function (complex) with the period r = 2/r/A^ 
defined by the desirable inter-channel spacing A^ • 

The first step of the example embodiment encoding algorithm deals with the specific 
representation of S(zy The sampling fimction is presented as the direct sum of LxN partial 
gratings, equally spaced in the fi:equency domain, and having relative phases (p^ and 

amplitudes dni—fdo, dj), n=l,...,N, l-J,...,L: 

Siz) = S i;^^ exp[i(2n + IN^l - 1) - 1 - NL)^ / 2 + )] . 

1-1 nsl 

The second step is to transform two-dimensional matrixes rf^, q>„i into a one- 
dimensional complex sequence m^y p = 1,2,..., JVL , where = d^^ expO>„/) . 
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Complete optimization of the sampling function would comprise numerical search for 
all jL X phases tp^ based on modified prior art techniques. For a high number of the recorded 
bits (e.g. LxN> 100) this approach may not be feasible due to excessive computational time. 
Therefoie» in the example embodiment the number of variables is decreased by presentmg the 
relative phases of the partial gratings in the form 9ni^<Pni » where for a fixed / are the 
N phases that optimize the /-th codeword. The phases ip^ are taken from a database of 
preliminary results, i.e. they are not calculated. 

Optixnization of all possible codewords (i.e. all possible combinations of {da, di) 
levels), for example using prior art methods, might be accomplished as a preliminary step for 
creation of the database. 

The rest of the phases, a set of L unknown variables (pf' , represents the relative phases 
between Ihe codewords. Below is the detailed description of how to obtain these L phases in an 
efficient way in the example embodiment. 

An algorithm to obtain L phases <pf has been developed, which comprises two distinct 
steps: minimization of the mean-square deviation of the sampling function, followed by an 
iterative procedure for further peak amplitude degradation. 

It has been shown that the peak amplitude of the sampling function is not high when its 
mean-square deviation is small. The mean-square deviation A(^) of the sampling function can 
be approximated accurately by an expression measuring the total energy of the aperiodic 
autocorrelation function Q ofsequrace m^: 



1 



^f2NL 



Elc.|" 



where 



= ' = l,2,...,Mr-l, 

and denotes complex conjugate. 
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The main advantage of this approximation is a high efiBciency of the nmnerical calculations 
when searching for the optimal set of phases. It is based on the possibility to derive an analytical 
expression for the gradient (direction towards the optimum) of the objective function, 



where Re stands for the real part of the expression. To calculate the above expression 
analytically we present the derivative of the autocoirelation function in the following form. 

We also present the autocorrelation index k and the sequence index p as 

k = ki + Nikj -1), p = Pi + NiPi -I) ' 

All possible values that indexes k md p can take in expressions (1) and (2) can be presented 
as three non-overlapping sub-sets: 

ifc, 6[l,L-l]U^i e[l,Ar|=>P2 "^jU/'i e[l + A„i\a 
ifcj e[l,L-l]U^^i e[l,Ar|=>P2 6[1 + *2,£]UPi ^[UNl 
K = L\Jky e [l,i\r-l] => Pa = A 6 [1 + ATI. 

For eadi of these sub-sets, calculation of the derivatives an^ ldq>j is a trivial exercise, i.e. 

(\ j ~ I 

Finally, the gradient of the mean-square deviation in the example embodiment can be 
calculated as 
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Im 



Ai+*i+^(y+*a-2) 



The above analytical expression forms a basis for the gradient method to find L phases 
qy'l , which minimize the mean-square deviation of sampling function Siz). One evaluation of 

both the mean-square deviation and its gradient requires in the order of N^L^ operations, i.e. 
exactly the same number of operations as in the case of uniform reflection spectrum (related 
art). However, a significant increase in the speed of tiie algorithm is caused by the decreased 
dimensionality of the phase space (fix>m LxN for related art methods toL for the method of the 
example embodiment). 

This has been found to produce optimised sampling functions characterized by low 
peaks in tiie refiractive index change. Further improvement is achieved in the example 
embodiment by using the second step of the optimisation procedure by a generalization of the 
iterative Gerchberg-Saxton algorithm or amplitude clipping. 

The Gerchberg-Saxton algorithm can be successfully used in cases when a small 
out-of-band response is not too crucial and the requirmients for ttie spectral resolution are not 
highly demanding. The key idea of the method is swapping between time/direction and 
firequency domains under the constraints that the amplitude of the complex sampling function is 
constant and the amplitudes of the central part of its spectrum are kept fixed to the desired 
levels. Iteratively, one translates virtually all modulations of the sampling function Siz) into its 
phase at the expense of appearance of small side-lobes in the spectrum of the envelope, the 
integral size of which is proportional to the mean-square deviation of 5'(z). 

An iterative clipping procedure is favourable in cases when an absence of side-channels 
in the reflection spectrum is essential. A complex error function is constructed by clipping S(z) 
at some level Si?. By subtracting the Fourier transform of the error function &om the finite 
snectrum of oricdnal Sfz^ and restoring the amplitude profile of the spectrum to the original 
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form (which includes setting the out-of-band response to zero), one decreases the maximum 
peak value of By gradually increasing level So^ one might significantly reduce the peak of 
the sampling function. 

The advantages of the above described Fourier space encoding approach include: the 
possibility to encode m-level data (m i.e. to write structures into a photosensitive medium 
which lead to the m-level reflection spectrum; and increased processing speed, with a potential 
for real-time optimization of high capacity storage devices (up to 10^ bits/mm^) based on the 
application of the developed block-coding approach. 

The example embodiment is extendable for bit sequences of arbitrary length. One can 
present a bit sequence of arbitrary length as NxL^xL^xL2^ ..X^^^ . In some cases, extra bits 
may be reqiiired to form ttie factorisation. Li an extension of the example embodiment, one first 
optimises NxL^ initial bits using the example embodiment described above. After that the 

algoritiun is repeated wifli x x £, bit sequence, but taking NxL^ as a "new" N and as 
a "new" . Then, one can keqp qtplying the algorithm, until the full bit sequence is optimised. 

■ 

In Figures 2(a) and (b), a sample two-level coding sequence in Fourier space for 675 
channels is shown to illustrate the potential of the present invention for storing a series of m- 
coded elements. The amplitude levels of the channels are 1.0 and 0.8. Figure 2(b) shows the 
detail of a portion of chaimels depicted in plot (a). Figures 2(c) and (d) depict the normalized 
(divided by V675) amplitude variation of the grating design presented on one period of the 
grating, corresponding to the bit structure of plot (a) based on the optimisation process 
embodying the present invention . Plot (d) depicts the clipped version of the of the data 
presented in plot (c). 

hi the following, an example embodiment of an optical data carrier writing system will 
be described. In Figure 3, an interferometric grating writing apparatus 10, under the control of a 
system controller 12, is utilised to write a sequence of transverse gratings (one- or multi- 
dimensional gratings) into a longitudinal strip-like waveguide. In the example embodiment a 
photosensitive optical strip 14 is continuously moved across an interference region 15 of the 
writing q[>paratus 10, as indicated by arrow 18. 

The writing apparatus 10, under the control of the controller 12, is arranged such that 
each grating is written transversely across into the optical strip 14. The strip may be continually 
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drawn past the writing beam in the writing process, in which case the writing apparatus 10 will 
have to be moved accordingly to write the gratings across the width of the strip 14. 
Alternatively the movement of the strip past the writing apparatus can be indexed such the at 
movement of the writing q)paratus during writing is not required. 

The optical strip 14, after having passed through the interference region 16 is then 
spooled into a disk-like shape as indicated by numeral 22. The strip 14 has a continuous 
reading face 20 along its length that lies generally perpendicular to the side of the strip 14 
through which the grating structure was written, and through which the gratings (e.g. 16) are 
read. The spooling of the optical strip 14 is done so that in the finished disk 24 the continuous 
reading face 20 of the strip 14 is accessible by a reading beam.. That is, the reading face 20 is 
parallel to the plane of the disk 22. 

The spooled disk 22 is then subjected to further processing (indicated at arrow 26) for 
producing a substantially solid disk 24. In the example embodiment, a refractive index 
matching optical glue can be used to adhere neighbouring windings of the coiled optical strip to 
form the disk 24. In another embodiment, the strip is made of a polymer material which is fiilly 
cured only from one side, and the grating structure is written only in that portion, with the top of 
the initially cured side forming the reading face. After spooling the strip into the disk-sh^e, the 
di^-shaped spool 22 is then cured further to achieve a complete solidification of the spool into 
a disk 24. In other embodiments, a possibility to un-roU the disk, e.g. to re-record it, may be 
provided for. In one such an embodiment, mechanical strength can be given to the spool 22 to 
form the disk 24 by encircling flie periphery of the spool in external plastic ring, the spool 22 
may also be supported on a substrate to provide adequate strength for handling or use. 

Turning now to Figure 4, it is noted that in the completed disk 24, the reading face 20 
lies within the plane of the disk 24. Accordingly, the disk 24 is suitsOjle for data reading e.g. by 
reflection, similar to the reading conducted in conventionally CD/DVD-type disk readers. 
Reading may comprise illuminating the grating structure using a non-collimated beam, e.g. a 
converging beam as shown in figure IB, so as to illuminate the grating structure at a plxirality of 
incidence angles. The reflected beam, which will only reflect at obtain angles, can then be 
detected by one or more detectors. The angles at which reflection occurs will correspond to the 
data encoded in the grating structure. It will be appreciated that the aspect of the present 
invention described with reference to the example embodiment in Figures 3 and 4 above, is not 
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limited to the use with Fourier-based encoding techniques, but can be utilised for any other data 
recording/storage approaches. 

It will be appreciated by the person skilled in the art that numerous modifications and/or 
variations may be made to the present invention as shown in the specific embodiments without 
departing fix>m the spirit or scope of the invention as broadly described. The present 
embodiments are, therefore, to be considered in all respects to be illustrative and not restrictive. 

For example, the partial grating sampling Amotions and thus the overall sampling 
function may be one- ore multi-dimensional. Also, the grating structures need not necessarily 
be UV written. Rather, they may be created by other techniques, including e.g. "stampm^* by a 
master strip or induced by other non-optical means such as e.g. pressure, chemicals, etc. 

In some embodiments, to separate the gratings firom each other in flie spooled or rolled- 
up state to avoid/reduce cross-talk between channels, some parts of the material strip may be 
left blank or UV exposed uniformly, i.e. those areas do not have any gratings written/induced. 

In the claims that follow and in the smnmaiy of the invention, except where the context 
requires otiberwise due to express language or necessary implication the word ^^comprising*' is 
used in the sense of 'including**, i.e. the features specified may be associated with fiirther 
features in varioiis embodiments of the invention. 
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